Resonant invisibility with finite range interacting fermions 
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We study the eigenstates of two opposite spin fermions on a one-dimensional lattice with finite 
range interaction. The eigenstates are projected onto the set of Fock eigenstates of the noninteracting 
case. We find antiresonances for symmetric eigenstates, which eliminate the interaction between two 
symmetric Fock states when satisfying a corresponding selection rule. 

PACS numbers: 34.30.+h, 05.30.Jp, 03.75.Lm, 05.45.Mt 



Introduction. Up to now a lot of interest in experi- 
ments on cold atoms has been focused on matter wave 
properties of the condensates which are described by the 
Hartree-Fock-Bogoliubov mean field model for weakly in- 
teracting quantum gases [1-5]. At the same time, the 
use of collision processes turns out to be a promising ap- 
proach to implement quantum gate operations [6[. A 
standard method for the description of such systems 
is to map them to Hubbard like lattice models where 
atomic physics provides a whole toolbox to engineer var- 
ious types of Hamiltonians for ID, 2D, and 3D Bose and 
Fermi systems. 

The interplay of interactions and discreteness leads to a 
set of interesting phenomena, including bound states, see 
e.g. [7UIH and |l2l-[l4|. In recent papers [H, [l6| we have 
studied properties of such bound states (also frequently 
coined quantum breathers) in various one dimensional 
Hubbard like models by considering two bosons or two 
fermions (with opposite spins) on lattices. The fermionic 
case adds to the complexity with the spin as an addi- 
tional degree of freedom. Consequently two fermions can 
form up to three different bound states, while two bosons 
form only one. In all these cases the interaction was as- 
sumed to be local, i.e. both particles interact only when 
occupying the same lattice site. In the present paper 
we consider fermionic particles with a finite range of in- 
teraction, as a more realistic description of experimental 
situations, which may be directly applicable in quantum 
computing, where the controlled interaction can be used 
to create entanglement with high fidelity. We analyze 
two particle eigenfunctions and identify resonance condi- 
tions for which two particles do not scatter despite the 
presence of a nonzero interaction. 

Model and spectrum. We consider one-dimensional 
periodic lattice with / sites described by an extended 
fermionic Bose-Hubbard (EFBH) model. The Hamilto- 
nian is given by 



H — Hq + Hjj + Hy , 
#0 = -^atviaj-i,* + 
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Here Ho describes the nearest-neighbor hopping, a =t,4- 
denotes the spin, Hjj and Hy describe the onsite and 
intersite (between adjacent sites) interaction between the 



particles with strengths [/, and V, respectively; 
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and 



are the fermionic creation and annihilation operators 
satisfying the anticommutation relations: [a,j i(T , a~^ a ,] = 

Sj,iSa,a', and [a ji(T ,ai i(T >] = [&+,., a^,] = 0. The sign of U 
and V is not specified. The Hamiltonian (pQ) commutes 
with the number operator TV = ^ hj whose eigenvalues 
are n = + n^, which is the total number of fermions 
in the lattice. In this work we focus on the simplest 
nontrivial case of n = 2, with = 1 and = 1. 

To describe quantum states, we use a number state 
basis |$ n ) = \ri1\n2" -rif) [l2[, where rti = + n^i 
represents the number of fermions at the i-th site of the 
lattice. The fermionic two particle states are generated 
from the vacuum \0) by successively creating a particle 
with spin down and spin up. 

The Hamiltonian (pQ) commutes with the translational 
operator T, which shifts all lattice indices by one. Its 
eigenvalues are r = exp(ik) with the Bloch wave number 
k = ^p, and 1/ = 0,1,2, ••• J -1 . 
Single- fermion states. In this simplest case, only one 
fermion is in the lattice (either with spin up or down) 



(n = 1), and the state is represented by \j) 
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\0). 



The interaction terms {Hjj and Hy ) have no contri- 
bution for a single particle. Thus the eigenstates of the 
Hamiltonian (pQ) are the eigenstates of Hq which are given 
by: 



The corresponding eigenenergies are 
£k = —2 cos(fc) . 



(5) 



(6) 



Two fermions. For the case of two opposite spin fermions 
(n = 2 with = 1 and = 1), each eigenstate is formed 



as a linear combination of number states with fixed n. 
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(7) 



For two particles, this involves N s = f 2 basis states, 
| $2)5 which is the number of ways one can distribute two 
fermions with opposite spins over the / sites with possible 
double occupancy. Then we define the basis state with a 
given value of the wave number k as in Ref. [16| and a 
complete wave function is: 

(/+i)/2 (/+i)/2 
l*2> = cil*i>+ E <*+l*i,+> + E ^-l*ii">- 

(8) 

Any vector in the corresponding Hilbert space is spanned 
by the numbers |ci, C2,+, C2,-, C3 5+ , C3 ? _ • • • ) and the vec- 
tors |3>i), and in two fermion case are de- 
fined as follows: 

l*0 = 77E(7)'"V,Al°>; 



i*i,+> 
1%-) 



77E(7)'" a MJ°>; 



(9) 
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We diagonalize the Hamiltonian (pQ) in the framework 
of the basis defined in (J5]) and derive the eigenenergies for 
each given Bloch wave number k from Hl^) = ^l^)- 
This leads to an / x / matrix whose elements H^j = 
2, ...,(/ + l)/2) are derived from 



= ^ = (*i,±|ff|$i) , H id = (* ii± |fr|$ j)± > . 

(10) 

We show in Fig. [T]the energy spectrum of the Hamilto- 
nian matrix ([TP]) obtained by numerical diagonalization 
for the case of opposite signs of interaction parameters 
U = 2 and V = — 3 and the form of the spectrum is 
similar to the one in Ref. [l6j|. Besides a two parti- 
cle continuum, three bound state bands are found. The 
eigenstates \&ki,k 2 ) of the continuum correspond to two 
fermions independently moving along the lattice as with 
zero interaction, and are derived from (|SJ). Their eigenen- 
ergies are given by : 



-4cos(fc/2) • cos(fci), 



(11) 



with k being the Bloch wave number and k\ = 2ttu/ (f — 
1) , being the canonically conjugated momentum of the 
relative coordinate (distance) between both particles and 
v = 0, ...,(/ — l)/2. Equation ([TT]) is the result of the 
sum of Bloch bands E± — — 2 cos(| d= k\) of two asymp- 
totically free particles [17j. 
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FIG. 1: Energy spectrum of two fermions of the EFBH chain 
with periodic boundary conditions for U = 2, V = — 3 and 
/ = 101. The lines follow from numerical diagonalization 
of the matrix [10] and symbols are the results of analytical 
computations for the bound states similar to the calculations 
in 



Weight functions in normal mode space. We transform 
to the basis of the symmetric and antisymmetric states 



V2 



V2 



(12) 

where a and s refer to the antisymmetric and the sym- 
metric states, respectively, j = 2, ...,(/ + l)/2. Note 
that |$i) is also a symmetric state. In this basis the 
matrix (|T0)) decomposes into two irreducible parts given 

by 



H s (i,j) = 



( U qy/2 
q*V2 V 
q* 



\ 



q* q 
q* p j 



(13) 



and 



(V q 

4* q 

q 



\ 



q 

q* -V ) 



(14) 



with q = 1 + r and p = r -(/+i)/2 + t -(/-i)/2 b xhe rank 
of the symmetric matrix is (/ + l)/2 and the rank of the 
antisymmetric matrix is (/ — l)/2. 

Our strategy is to compute an eigenstate for the inter- 
acting case, and use a weight function to expand it in the 
basis of the eigenstates of the noninteracting case. For 
this purpose we fix the Bloch momentum fc, and choose a 
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condition for zero weight: 



FIG. 2: Weight function for two seed mode numbers ki = 
0.237T and k\ = 0.797T, with the onsite and intersite interaction 
parameters U = 0.1 and V — 0.08. Here k — 0.5tv and / — 
101. Dashed lines are the results using formula (|16p 



seed eigenstate ) of the unperturbed case with seed 

mode number k. Upon switching on the interaction it 
becomes a new eigenstate |^ r ^ 1 ), which will have over- 
lap with several eigenstates of the unperturbed case. We 
expand the eigenfunction of the perturbed system using 
first order perturbation approximation: 



1 k\ 



ki 



-pT-^^kO- ( 15 ) 



From expansion (Tl5| it follows that the off-diagonal (fci 7^ 
k\) weight function at the first order is given by : 



C s (fci;fci) = 



{2U\q? + VElEl?\{^l i ) 



4|g|« 



(16) 



with \q\ 2 = 2 + 2cos(/c) and and E~_ are the eigenen- 

ergies of the unperturbed system given by Eq. ([IT]) . 

Symmetric states. First we consider the case of small 
interaction constants and as expected we find localization 
in normal mode space. For instance, in the case of dom- 
inant onsite interaction constant U = 1 and V — 0.1, 
quite similar results to those obtained in Ref. [15| are 
derived and the result of the perturbation formula (fTBj) 
matches pretty well with those of the diagonalization pro- 
cedure. 

But if now we take interaction constants U, V with 
comparable values (but again in perturbative limit) 
one will unavoidable deal with additional antiresonance 
structure presented in Fig. [2j where the weight function 
vanishes exactly. The appearance of these structures fol- 
lows from the analytical formula (fT6|) . Indeed, in per- 
turbative limit of interaction constants one can find such 
a seed k\ and probe k\ wavenumbers that the weight 
function becomes exactly zero. We find the following 



U 
V 



-2 cos ( &i) cos (&i 



(17) 



It also follows from eq. ([TT)) that there is a critical wave 
number given by 



k x = arcos{- — ) . 



(18) 



An antiresonance appears only if the following inequali- 
ties are satisfied: n — kl<ki<kl (here for simplicity we 
assume both interaction constants positive). As it is seen 
from Fig. [2] the perturbative limit ([T6]) works well even 
in case of presence of an antiresonance. Equation ([T7]) 
further tells that an antiresonance will be observed even 
for U = 0. In this case the seed wavenumber k\ = tt/2 
is not modified by interaction V. On the other hand if 
V — the antiresonances are not observable. 

For larger values of interaction constants the perturba- 
tive predictions will get significant corrections. To show 
this we plot three dimensional graphs of weight function 
versus seed k\ and probe k\ wave numbers for various 
values of interaction constants in Fig. [3j As seen for 
small values of interaction constants the track of the an- 
tiresonances keeps the symmetry in the seed-probe mode 
number space traces predicted by perturbation theory. 
However for large interaction constants this symmetry is 
lost. 

Antisymmetric states. The structure of the antisym- 
metric matrix ([T4]) suggests that the weight function for 
antisymmetric states can be computed as: 



C a (krM) 



K~ E V 2 



(19) 



and according to this formula the weight function does 
not develop antiresonances. This has been confirmed by 
numerical diagonalization. 

Discussions. Let us discuss the meaning of the ob- 
served antiresonances. Two particles, when travelling 
along the lattice, will meet, interact, and scatter. If pre- 
pared in an initial symmetric noninteracting seed state, 
the particles will scatter into all other available noninter- 
acting symmetric states - except for one special. This is 
because the scattering can go either via the onsite inter- 
action U or via the intersite interaction V. A correspond- 
ing destructive interference makes the amplitude in this 
particular scattering state exactly zero. Antisymmetric 
states have strict zero occupation on the same site, and 
therefore only one scattering path (using V) is left. Con- 
sequently they do not show antiresonances. But they 
will, if we add even more distant (e.g. next-to-nearest- 
neighbor) interactions. 
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